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Abstract 

We calculated the dynamics of nuclear wave packets in coupled electron-vibration sys- 
tems and their nonlinear optical responses. We found that the quantized nature of the 
vibrational modes is observed in pump-probe spectra particularly in weakly interacting 
electron-vibration systems such as cyanine dye molecules. Calculated results based on 
a harmonic potential model and molecular orbital calculations are compared with ex- 
perimental results, and we also found that the materials parameters regarding with the 
geometrical structure of potential energy surfaces are directly determined by accurate 
measurement of time-resolved spectra. 

1 Introduction 

Optical control of quantum-mechanical states has attracted attention of many authors[lJ. In 
particular, recent progress in laser technology has made it possible to synthesize desired optical 
pulses, and they are expected to be a tool to control physical properties of various materials, 
e.g., luminescence intensity, lifetime of excited states, or chemical reaction rate. 

We, however, have not found effective methods for coherent control of those properties, since 
it is impossible to design appropriate optical processes without the knowledge of the structure 
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of the electronic states and/or the potential energy surfaces (PESs) of materials. Thus, it is 
important to determine the values of the material parameters regarding with the ground-state 
and/or excited-state PESs to realize coherent control of above-mentioned physical properties. 
As a first step toward this approach, we study methods to trace the motion of wave packets on 
each PES after photoexcitation in order to understand the structure of PESs. 

Experimental studies on wave packet dynamics have been performed with ultrashort laser 
pulses since the pioneering work by Zewail et al.[2\, and previous studies on molecules in solution 
give a clue to understand the motion of wave packets created by optical pulses[3l HJ EE El- 
Alt hough the nonlinear response of wave packets on PESs has been understood by means of 
the motion of a "classical point" on PESs, our recent experiments on cyanine dye molecules 
have shown that the wave packet dynamics in those molecules is understood only by taking into 
account the quantum nature of the vibrational modes of the molecules|7j. According to our 
preliminary calculations [7J, this particular feature is obtained when the interaction between 
electrons and the vibrational modes is weak. To be more precise, the quantum nature of 
vibrational modes becomes prominent when the average number of vibrational quanta contained 
in each wave packet is ~ 1, which shows that the wave packet motion is restricted in the vicinity 
of the bottom of the PESs. In this case, the relaxation energy of the Franck-Condon state is 
comparable to the energy of vibrational quanta hu>. 

In this paper we calculate the dynamics and the nonlinear optical response of wave packets 
on excited-state PES when quantization of vibration modes plays an important role, i.e., in 
weakly coupled electron-vibration systems. In this particular case, we found that the impulse 
approximation is effective to analyze the experimental data. We also show that accurately 
measured pump-probe signals help us determine the values of the material parameters regarding 
with the structure of PESs, which are, as mentioned before, important to design the coherent 
control methods of wave packets by laser pulses. 

The paper is organized as follows: in Section [2] we introduce the model and its basic prop- 
erties. We show the calculated results of the pump-probe signal including the molecular orbital 
calculation for cyanine dye molecules in Section [3J Section H] is devoted to discussion and 
conclusions. 

2 Model 

In order to describe the dynamics of wave packets in typical organic molecules, we take into 
account relevant two electronic states which correspond to the ground and the excited states 
and interact with multiple vibrational modes. Thus, we employ the simplest model suitable for 
the present purpose, i.e., a two-level model on potential surfaces with M modes described by 
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the following Hamiltonian: 




In Eq. ([!]) | T)) denotes the excited electronic state of a single molecule, and its ground electronic 
state is | j)). Pj an d Qj are the normalized momentum and the coordinate of the vibration modes 
of the molecule, and Sj is the dimensionless coupling constant between an electron and the j-th 
vibrational mode. 

Since, in weakly coupled electron-vibration systems, the wave packet motion is restricted 
only in the vicinity of the bottom of PES as mentioned before, the parabolic approximation is 
appropriate in the present study In particular, a Gaussian wave packet which is generated by 
the vertical transition by optical field does not distort while it moves on parabolic PESs. This 
property appears to be quite useful, as the optical response within the impulse approximation 
is obtained analytically as shown in the next section. 

In the rest of the paper, we assume that the vibrational frequencies on both PESs are the 
same, i.e., we take tj = 0, which is appropriate for cyanine dye molecules as shown in the 
next section. The details of the quantization of the Hamiltonian ([1]) is also described in the 
Appendix El and we show in this section only the eigenvalues and eigenstates of Hamiltonian 
([T]) corresponding to the ground-state PES and the excited-state PES for tj =0: 
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where the boson number states (Fock states) are defined by 
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3 Calculated results 

3.1 Pump-probe signals of wave packets on harmonic potential sur- 
face 

When we consider that the temporal width of both the pump and the probe pulses are small 
compared with the period of vibration modes 2n/ujj, we apply the impulse approximation to 
calculate the probe signal [8] as 
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where the time-dependent polarization induced by the probe pulse P(T,t) is given by 



P(T,t) = l -{®{T)\[V{t - T), V]\HT))0(t - T). (8) 

Q and T denote the frequency and the delay of the probe pulse, respectively, and 6{t) is the 
step function. 

The wavefunction \$>(T)) is the Franck-Condon state multiplied by the time-evolution op- 
erator with regard to TC and is described by 

|$(T)> = e~— |si,s 2 ,...,sa/;T) 

= e- ieT \ Sl e- iuJlT , s 2 e-^ 2T , ...s A/e -^ r ; |), 

(9) 

where |oti, a 2 , ■■■otM'i T) is a Gaussian wave packet on the excited-state PES, i.e., a coherent 
state of the vibration modes defined by 
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The polarization operator is given by V = | T))(U I + I l))((T I an d its time-dependent form 
V(t) is 

V(t) = e^Ve 

in the interaction representation. 

After algebraic calculations, P(T,t) is obtained as 
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P(T,t) = Im 
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Combining Eqs. ((7j), ({12]) . and f|T4l . we found that finite probe signal is obtained only when 
the resonance condition 
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is satisfied. 

For M = 1, P(i, T) has a far simpler form and A(Q; T) is analytically obtained for arbitrary 
value of Q and T. When we consider the case for M = 1 the indices of the parameters are 
omitted for brevity in the rest of the paper. The Fourier transformation of P(t, T) is performed 
in a straightforward manner and A(Q; T) is given by 
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where J n (z) is the Bessel function of the first kind and [a] denotes the largest integer which 
does not exceed the value of a. When s 2 ^> 1, vibration eigenstates are densely distributed 
between the Franck-Condon state and the bottom of the excited-state PES, which means that 



5 




n=5o 



n=i5o 



a=250 



n=350 



(b) s=10 



£2=450 



£1=550 



delay(l/a>) 



5 

delay(l/to) 



10 



Figure 1: Probe response for M — 1 and (a) s — l,e — 10, (b) s — 10, e 
probe frequency in units of hu. 



100. Q denotes the 



the classical picture based on continuous PESs is recovered. For example, A(Q; 2nn / u) (n = 
0, 1, 2, 3...) is obtained as 

A(n;2nir/uj) oc .? 5{tt-e-ju). (17) 

The coefficients on each 5-function in Eq. f|T7|) show that the probe response as a function of 
probe frequency is given by the Poisson distribution which approaches the normal distribution 
in the limit of s — > oo. 

We numerically calculated A(Q;T) for arbitrary values of Q or T, and Figs. [T]-(a) and (b) 
show the probe response as a function of T for several values of O for s — 1 and s = 10, 
respectively. When s = 10, the temporal behavior of the probe response is understood by 
means of a classical PES, i.e., finite reponse to the probe pulse is observed only when the 
center of the wave packet passes the point on the PES at which the probe pulse resonates to 
the transition between | j) and | |). However, when s — 1, the quantum nature of the vibration 
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mode appears to be important and the shape of the probe spectra as a function of pump-probe 
delay T apparently differs from that obtained from classical PESs, i.e., the probe response is 
finite even when the center of the wave packet is not on the resonating point of the PES to the 
probe pulse. In other words, the probe signal is temporally broadened due to the finite spatial 
width of the quantum- mechanical wave packet. 

Reflecting the motion of the wave packet, the probe signal shows a oscillatory behavior as 
a function of T. In the classical limit the peaks of the oscillation in the spectra gradually move 
as the probe frequency is changed, whilst the peaks and the valleys in the spectra suddenly 
interchanges with each other as the change of the probe frequency in the quantum-mechanical 
regime. Thus the time-resolved spectra translationally move along the temporal axis as if their 
phase shifts by it. As Eq. (j!4p shows, however, the phase of the oscillation as a function of 
T does not alter, and thus the peak shift in the spectra cannot be discussed in terms of the 
phase of the oscillation. Hence, we stress that time-dependent shift of the phase in each Fourier 
component which is often taken into account in experimental studies is not necessary to analyze 
the time-resolved spectra to discuss the wave packet motion. 

The quantum nature of the PESs is also shown in observing the "turning points" of the wave 
packet motion. When probe frequency is varied, we observe spectral structure which shows the 
turning points and they help us determine the values of material parameters [7J. Since, however, 
both the average and the variance of the number of vibrational quanta in a Gaussian wave packet 
\s t) are s 2 , the probe response of wave packets is obtained for higher/lower probe frequencies 
than those at the classical turning points. As a result, the spectra become blurred as shown 
in experimentsjTj. However, the peaks of the spectra still gives the turning points in classical 
sense, and therefore we expect that accurate measurement of time-resolved spectra is still an 
effective tool to obtain the values of materials parameters. 

It has been pointed out that the probe signals include the optical response of wave packets 
on the ground-state PES as well as those on the excited-state PES [I]. Since, as shown in Eq. 
(J7J), the impulse approximation takes into account the signals only from the excited-state wave 
packets, it is necessary to consider the validity of the impulse approximation for the probe 
pulses with finite width. 

Using the standard third-order perturbation theory [8], we found that the contribution of 
ground-state wave packets to the probe response can be estimated by the transition probability 
between the ground and excited PESs, 

# r 2, ^{2S?(1- COS WiT)}"* . . 

r nm 2 ...n M = 11 exp{-2s j (l - coswj-r)} J - j , (18) 

j=i n r 

where r denotes the temporal width of the pump pulse. The detail of the calculation to derive 
Eq. (fl8|) is presented in Appendix [Bl 
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Figure 2: Transition probability r n (n = 0, 1, 2, 3, 4) for M — 1 and (a)s = 1, (b)s = 10. 

Figure [2] show r n (n = 0, 1, 2, 3, 4) for s = 1 and s = 10 for M = 1. For s = 1, r n for n^0 
is small compared with r for out < 0.5, which corresponds to the case where the pump pulse 
is shorter than the tenth of the vibrational period. For this short pulse, the ground-state wave 
packets generated by the pump pulse are dominated by |0 j) which does not contribute to 
the temporal oscillation of the probe spectra, and hence only the motion of excited-state wave 
packets is reflected on the time-resolved spectra. 

On the other hand, for s = 10, r n for n ^ is comparable to tq and hence the ground- 
state wave packets consist of multiple vibrational states. In this case, the time-resolved spectra 
contain oscillatory structure due to the motion of ground-state wave packets as well as that of 
the excited-state wave packets, and the impulse approximation fails. 

Hence, we found that the validity of the impulse approximation is estimated by calculating 
the values of r n . Comparing r with r\, we found that the impulse approximation works when 
p c = ri/r = 2s 2 (l — cos cur) or, in general, all of the corresponding factors pj, for j-th mode 
are sufficiently small. In particular, for weakly coupled electron-vibration systems such as dye 
molecules (s < 1), the probe signals given by the wave packets on the ground-state PES are 
independent of the delay T for short probe pulses used in experiments, and thus the present 
calculation is valid to discuss the motion of the excited-state wave packets. 

3.2 Pump-probe response of cyanine dye molecules 
3.2.1 Molecular Orbital Calculations 

Before we apply the present theory to the optical properties of DTTCI (diethylthiatricarbocya- 
nine iodide) measured in our previous study [7], we calculate the electronic state of DTTCI by 
the molecular orbital calculations. DTTCI is dissociated into C25H25N 2 S^ and I~ in the ethanol 
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solution. Therefore, we consider DTTC + (C 25 H25N 2 S2~) as a model molecule, where the solvent 
effect is not considered. The optimized geometry for the singlet ground state of DTTC + was 
obtained in the restricted Hartree-Fock (RHF) approximation, and the vibrational frequency 
analysis was performed. The optimized geometry for the first singlet excited state was ob- 
tained in the singly excited configuration interaction (CIS) approximation, and the vibrational 
frequency analysis was performed. The basis set employed in this work was 6-31G(d)[9j, and 
all calculations were performed with the Gaussian 98 program package [9J. 

The probability amplitudes for the highest occupied molecular orbital (HOMO) and the 
lowest unoccupied molecular orbital (LUMO) of DTTC + are given in Fig. [31 The character of 



Figure 3: The probability amplitudes for HOMO and LUMO of DTTC+. 

the first excited state in the CIS calculation is the excitation from HOMO to LUMO. The large 
oscillator strength for this excitation, 2.94 in the CIS calculation, indicates that DTTCI is the 
efficient dye molecule. 

Schematic representation of the potential energy surfaces for the ground (E g (R)) and the 
first excited states (E X (R)) of DTTC + is illustrated in Fig. HI The minimum energy configu- 
ration for the ground state is shown as A, and that for the first excited state is shown as B. 
The difference of the configuration between A and B is quite small, and the energy difference 
E X (A) — E X (B) is so small as 52.5 meV. The difference between the vibrational frequency for 
the first excited state and that for the ground state is also quite small. Therefore, the weak 
coupling model can be applied to this problem, and the vibrational modes for the first excited 
state can be considered to be the same as those for the ground state. 

In the Franck-Condon picture, absorbing the pump light, the DTTC + molecule is vertically 
excited from the ground state to the first excited state, and the DTTC + molecule is placed 
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Figure 4: Schematic representation of the potential energy surfaces for the ground and the first 
excited states of DTTC + . 



on the configuration A on the potential energy surface of the first excited state (E X (A)). 
In order to investigate the time-dependent nuclear dynamics of DTTC + molecule after the 
excitation, molecular dynamics (MD) calculation is necessary. Considering the nucleus as 
classical particles, the MD of DTTC + can be described by the Newtonian equation of motion 
with the potential (^(R)). Although we have not performed the full MD calculation, we have 
calculated the force at (^(A)), 

dE r (R, 

(19) 

R=A 
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as the measure of the initial dynamics. The vector F has 162 components corresponding 
to the x, y, and z directions of each atom. Then, the scalar product of F with each unit 
displacement vector Ui(i = 1 ~ 156) corresponding to each harmonic vibration on (E X (R)) has 
been calculated. The vector u« has 162 components as well. 

The vibrational frequency modes, having larger absolute value for the scalar product with 
F, are given in Table [IJ where F is normalized. The mode 45 has the largest absolute value for 
the scalar product. Considering the beat of our experiment has the period of ~ 230 fsec, the 
frequencies for mode 15 is close to the experimental frequency. For further verification of the 
experiment, we have to perform the MD calculation . 



3.2.2 Pump-probe spectra of DTTCI 

To compare with the experimentally observed pump-probe spectra[7], we calculated the probe 
response based on Eq. ([121 for M — 2. The values of the parameters of the vibrational 
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0.11 


130 


1785.3942 


0.08 


131 


1827.9013 
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Table 1: The vibrational frequency modes on Ex(R), having larger absolute value for the scalar 
product with F. 
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frequencies are hu>i = 18meV (mode 1) and hu)2 = 58meV (mode 2). The electron- vibration 
coupling constant for the first mode s% is fixed to 1. These values are taken from Ref. [7]. 

The calculated results of the pump-probe spectra for S2 = 0.6 is shown in Figs. 13- (a) and 
(b). Although probe response is finite only for discrete values of probe frequencies (Eq. [T5l) . we 
draw lines which connect spectral lines with finite response, so as to make the overall behavior 
of spectra clearly shown. 

The motion of the excited-state wave packet is understood by tracing the peak of the spectra. 
However, as shown in Fig. [3(b), the trajectory of the peak of the spectra corresponds to that for 
the vibration mode 1, although fine structure due to the mode 2 is also seen in the spectra (see 
Fig. [3 (a)). In fact, the amplitude of the oscillation in Fig. 13(b) is ~ 70meV which corresponds 
to the distance between the two turning points for the mode 1 (Asfhui = 72meV). 

As the value of Sj decreases, the trace of the peak of the spectra is blurred and the turning 
points become unclear as we mentioned in Section 13.11 and for this reason the turning points 
corresponding to the mode 2 are not seen in the spectra. As Fig. [3(a) shows, however, the 
effect of the mode 2 is clearly seen in the fine structure of the spectra, which means that we 
understand the wave packet motion corresponding to the mode 2 by accurate measurement of 
time-resolved spectra. To be more precise, we found that the shape of the spectra is sensitive 
to the value of S2 in the present case and that the S2 is determined by careful observation of 
the spectral shape. In Figs. [6]-(a)-(c) we show the results for s 2 = 0.5, 0.6, and 0.7 with the 
values of the other parameters fixed to those in Fig. [3 The figures show that the shape of 
the spectra changes rapidly as s 2 is varied. Thus, we can determine the materials parameters 
regarding with weakly coupled mode by a sort of pattern matching methods for time-resolved 
spectra, though the precise methods to determine these values have not been well-established. 
In the case of DTTCI, we found that s 2 = 0.6 is most feasible by using Figs. [3(a)-(c). Hence, 
the relaxation energy of the excited wave packet is 39meV, which is consistent with the results 
of molecular orbital calculation shown in Section 13.2.11 

4 Discussion and conclusions 

In this paper we calculated the nonlinear optical response of wave packets on quantized potential 
surfaces. We found that, due to the discreteness of the vibrational eigenstates, the probe signals 
as a function of pump-probe delay suddenly changes its behavior following the stepwise change 
of the energy levels of the molecule. This aspect is particularly important when electron- 
vibration interaction is weak, i.e., the wave packet motion takes place in the vicinity of the 
bottom of the excited-state PES. In this case, a simple parabolic approximation for PESs is 
quite good and the analytical formula of the probe signal within impulse approximation is 
sufficient to discuss experimental results. 
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Figure 5: Calculated pump-probe spectra for DTTCI as a function of probe photon energy and 
the pump-probe delay time: (a) bird's-eye view (b) palette-mapped plot. The values of the 
parameters are: Twj\ = 18meV, hu>2 = 58meV, he = 1.5eV, Si = 1, and s 2 = 0.6. 
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Figure 6: Calculated probe response at Q = 1.52eV for huJi = 18meV, T1UJ2 = 58meV, Si = 1 
and (a)s 2 = 0.5, (b)s 2 = 0.6, (c)s 2 = 0.7. 
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The calculated results agree well with the pump-probe experimental results on DTTCI, and 
it is shown that we can estimate the strength of electron-vibration interaction directly from 
experimental results. These results are also consistent with the molecular orbital calculations. 
Thus, we stress that the accurate measurement of pump-probe signals combined with theoretical 
calculations is a powerful tool to understand the wave packet dynamics clearly. Once the 
dynamics of the wave packets are sufficiently known, we try to design the control method of 
them, which will bring us, e.g., the efficient way for chemical reaction control, optimization of 
luminescence intensity from biological samples, or accurate control of qubits stored in various 
materials such as quantum dots or molecules. We mention that, as a typical example of the 
control method, the initial condition for the excited wave packet is controlled by pulse chirp [10, 
ITTl [T2"] and that the probe response is modulated by the presence of chirp. We, however, 
also note that those results have not reached to control the wave packet motion observing the 
quantum-mechanical wave packets and hence the design of the control in coherent regime is 
still on the way even in simple molecular systems. 

We should mention that the present calculation does not take into account the effect of 
decoherence of the vibrational states. Since we are interested in molecules in solution, the 
relaxation of excited wave packets is important to understand their physical properties even 
in sub-picosecond regime. These effects necessarily involve the mixed-state description of the 
wave packets and coarse graining of the dynamical variables consisting of the "environment" . 
It is necessary to obtain various knowledge of the total system (molecules+environment) to 
develop a quantitative theory of the decoherence in the present system, and thus the problem 
of decoherence is left for the future. We, however, stress that the present results are still useful 
as a basis of the study of wave packet dynamics even when the relaxation processes are of 
importance, 
acknowledgements 
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A Quantized form and the eigenstates of the Hamilto- 
nian 

After canonical quantization of the Hamiltonian (jTJ, it is rewritten as 
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where the creation operator of the vibration modes of the ground-state PES is 

t _ F„, s p i 



It is diagonalized by introducing unitary transformation 

M 

U = Y[Ui, (22) 
i=i 

where 

Ui = exp |^^((ai) 2 - a 2 )} exp {^"(4 - a*)} • (23) 

The eigenvalue and the corresponding eigenvector of the Hamiltonian on the excited-state PES 
are given by 
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and the creation operator of the phonon is 
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where 

e «* = (1-^-1/4 (27) 

When the frequency of the vibrational modes for two PESs is the same, i.e., ti = 0, the above 
results are identical to those previously obtained [T3]. 

The eigenstates on the excited-state PES for each mode U}\rii) are expressed by 
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Thus, the matrix elements of the dipole operator V is given by 
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{(ajpaf} denotes the symmetrically ordered product of the operator in the brace[I3], and the 
sum over k is taken only for the integral value of mi+ 2 ~ 3 . 



B probe response by ground-state wave packets 

Amongst the third-order polarization terms which appear in the standard perturbation theory [8], 
the contribution of ground-state wave packets to pump-probe signals is described by, for exam- 
ple, 

fj 3) (*) = f f f dt 1 dt 2 dhE 1 {t l )E l {t-t l -t 2 )E 2 {t-t 1 -t 2 -h){Q | \VV{ti)V{ti+t 2 )V(ti+t 2 +t 3 )\0 I), 
M Jo Jo Jo 

(31) 

where Ei(t) and E 2 {t) denote the electric field of the pump and the probe pulse, respectively. 

When the pump and the probe pulse are short and they are separated on the temporal axis 
sufficiently, Eq. (13B for M — 1 becomes 

if) ( t ) „ §lS 2 £ e m " T (se~^; | \V\n |)(n | iVe'^^VlO |), (32) 

n=0 

where Si and S 2 are the pulse area of the pump and the probe pulses, respectively, r and T 

i7i(t—T) 

denote the duration of the pump pulse and the pump-probe delay. Since (0 j \Ve s V\0 j) 
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corresponds to the ground state absorption which is independent of T, the vibronic ground 
state |0 |) in Eq. (I3"2"j) does not contribute to the T-dependent part of Pjj 3 \t), and is irrelevant 
to the oscillatory behavior of the time-resolved spectra. Hence, the magnitude of the temporal 
oscillation in time-resolved spectra given by the ground-state wave packets is estimated by the 
factor | (se~ llJJT ] f |V"|7& j) | 2 which is interpreted as the transition probability of vertical transition 
between excited-state wave packets and ground-state wave packets. 

It is quite easy to extend the above discussion to mult i- dimensional cases (M 7^ 1), and we 
obtain the generalized formula to estimate the effect of ground-state wave packets as 

\( ni n 2 ...n M I \V\ Sl e-^ T s 2 e-^ T ...s M e-^ T -V\ 2 

¥t r i, „ {2sf(l -coster)}? . , 

TT exp{-2s?(l - coauir)}^-^ ; — 3 —^-. 33 



T n\ri2...riM 
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